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COBE-DMR-Normalized Dark Energy Cosmogony 

Pia Mukherjee 1,2 , A. J. Banday 3 , Alain Riazuelo 4 , Krzysztof M. Gorski 5,6 , and Bharat Ratra 1 

ABSTRACT 

Likelihood analyses of the COBE-DMR sky maps are used to determine the normalization 
of the inverse-power-law-potential scalar field dark energy model. Predictions of the DMR- 
normalized model are compared to various observations to constrain the allowed range of model 
parameters. Although the derived constraints are restrictive, evolving dark energy density scalar 
field models remain an observationally-viable alternative to the constant cosmological constant 
model. 

Subject headings: cosmic microwave background — cosmology: observation — large-scale struc- 
ture of the universe 



1. Introduction 

Indications are that at the present epoch the universe is dominated by dark energy (see, e.g., Peebles & 
Ratra 2003 for a review and Bennett et al. 2003a for a summary of the recent WMAP results). Einstein's 
cosmological constant A is the earliest example of dark energy, and more recently scalar field models in 
which the energy density slowly decreases with time, and thus behaves like a time- variable A, have been the 
subject of much study (see, e.g., Peebles 1984; Peebles & Ratra 1988, 2003; Padmanabhan 2003). 7 

In this paper we focus on a simple dark energy scalar field (c/>) model. In this model the scalar field 
potential energy density V(<j>) at low redshift is oc <j)~ a , with a > (see, e.g., Peebles & Ratra 1988; Ratra & 
Peebles 1988). Consistent with the observational indications, the cosmological model is taken to be spatially 
flat (see, e.g., Podariu et al. 2001b; Durrer, Novosyadlyj, & Apunevych 2003; Page et al. 2003; Melchiorri & 
Odman 2003). 
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Dark energy models, which assume an early epoch of inflation to generate the needed initial energy 
density fluctuations, do not predict the amplitude of these fluctuations. At present, the most robust method 
to fix this amplitude, and hence the normalization of the model, is to compare model predictions of the large 
angular scale spatial anisotropy in the cosmic microwave background (CMB) radiation to what is observed. 
To this end we compute the model predictions as a function of the parameter a and other cosmological 
parameters by following Brax, Martin, & Riazuelo (2000). We then determine the normalization amplitude 
by comparing these predictions to the COBE-DMR CMB anisotropy measurements (Bennett et al. 1996) 
by following Gorski et al. (1998). 

In linear perturbation theory, a scalar field is mathematically equivalent to a fluid with time-dependent 
equation of state parameter w = p/ p and speed of sound squared <? s = p/p, where p is the pressure, p the en- 
ergy density, and the dots denote time derivatives (see, e.g., Ratra 1991). The XCDM parametrization of this 
dark energy model approximates w as a constant, which is accurate in the radiation and matter dominated 
epochs but not in the current, dark energy scalar field dominated epoch. This XCDM approximation leads 
to particularly inaccurate predictions for the large angular scale CMB anisotropy. We emphasize, however, 
that we do not work in the XCDM approximation, rather we explicitly integrate the dark energy scalar field 
and other equations of motion (Brax et al. 2000). 

With a robust determination of the normalization of the inverse-power-law potential dark energy scalar 
field it is possible to test the model and constrain model parameter values by comparing other model 
predictions to observational data. In this paper we compare the predicted value, as a function of model 
parameters, of the rms linear mass fluctuation averaged over an 8 h~ x Mpc sphere 8 , <5M/M(8ft. -1 Mpc) or 
(Ts, to observational estimates of this quantity. We also present qualitative comparisons of the predicted 
fractional energy-density perturbation power spectrum P(k) to various observations. 

Such tests probe a combination of both the early universe physics of inflation and late universe dark 
energy physics. 9 Other neoclassical cosmological tests arc largely insensitive to the early universe physics of 
inflation and are hence more direct probes of dark energy. Of special interest are tests based on gravitational 
lensing (see, e.g., Ratra & Quillen 1992; Waga & Frieman 2000; Chae et al. 2002), Type la supernovae 
redshift-apparcnt magnitude (see, e.g., Podariu & Ratra 2000; Waga & Frieman 2000; Leibundgut 2002; 
Tonry et al. 2003), and rcdshift-angular size (see, e.g., Chen & Ratra 2003a; Podariu et al. 2003; Zhu & 
Fujimoto 2003, also see Daly & Djorgovski 2003) data. 

In § 2 we summarize the techniques we use to determine the DMR estimate of the CMB rms quadrupole 
moment anisotropy amplitude Qrms-PS for the inverse-power-law potential dark energy scalar field model. 
Results for Q rms -ps and large-scale structure statistics for the DMR-normalized models are given in § 3. 
These statistics are compared to current observational measurements in § 4. Our results are summarized in 
§5. 



8 Here the Hubble constant H = WOh km s _1 Mpc -1 . 

9 Another such test is the detection of the late integrated Sachs- Wolfe effect from the cross-correlation of various large scale 
structure maps with maps of CMB temperature anisotropy (see Boughn & Crittenden 2003; Fosalba, Gaztanaga, & Castander 
2003; Myers ct al. 2003; Nolta et al. 2003; Scranton et al. 2003). 
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2. Summary of Computation 

We consider a scale-invariant primordial energy-density fluctuation power spectrum (Harrison 1970; 
Peebles & Yu 1970; Zel'dovich 1972), as is generated by quantum fluctuations in a weakly coupled scalar 
field during an early epoch of inflation (see, e.g., Fischler, Ratra, & Susskind 1985), and consistent with 
the observational indications (see, e.g., Spergel et al. 2003). More precisely, the fractional energy-density 
perturbation power spectrum we consider is 

P(k) = AkT 2 {k), (1) 

where k (0 < k < oo) is the magnitude of the coordinate spatial wavenumbcr, T(k) is the transfer function, 
and A is the normalization amplitude defined in terms of the Bardeen potential. 

The CMB fractional temperature perturbation, ST/T, is expressed as a function of angular position, 
(6, 4>), on the sky via the spherical harmonic decomposition, 

„„ oo £ 

— (e,cj>)=J2 E a tm Y lm {e,4>). (2) 

£=2 m=-£ 

The CMB spatial anisotropy in a Gaussian model 10 can then be characterized by the angular perturbation 
spectrum Cg, defined in terms of the ensemble average, 

(ai m ae/ m i*) = Cedu'Smm', (3) 
where the S's are Kronecker delta functions. 

The C/s used here were computed using the Boltzmann transfer code of Brax et al. (2000). The 
computations here assume a standard recombination thermal history, and ignore the possibility of early 
reionization, tilt, gravity waves, or space curvature. While observed early reionization (Bennett et al. 2003a) 
does affect the Ce's on the large angular scales of interest here, the effect is not large and an accurate 
quantitative estimate awaits a better understanding of structure formation. At present there is not much 
observational motivation for the inclusion of tilt, gravity waves, or space curvature. 

The dark energy scalar field model we consider here is characterized by the potential energy density 
function 

V(<f>) = n<b- a , (4) 

where the constant k has dimensions of mass raised to the power a + 4 (Peebles & Ratra 1988). In the limit 
where the exponent a approaches zero, the scalar field energy density is equivalent to a constant cosmological 
constant A. Here the exponent a > 0. 

We evaluate the CMB anisotropy angular spectra for a range of inverse-power-law scalar field potential 
exponent a spanning the interval between and 8 in steps of unity, a range of the matter density parameter 
ilo spanning the interval between 0.1 and 1 in steps of 0.1 and for fio = 0.05, and for a variety of values 
of the hubble parameter h and the baryonic-mass density parameter Qb- The values of h were selected to 
cover the range of ages consistent with current requirements (to = H Gyr, 13 Gyr, or 17 Gyr, see, e.g., 



10 Simple inflation models are Gaussian, with fluctuations generated in a weakly coupled field (see, e.g., Ratra 1985; Fischler 
et al. 1985), consistent with observational indications (see, e.g., Mukhcrjcc, Hobson, & Lasenby 2000; Park et al. 2001; Komatsu 
et al. 2002, 2003; Santos et al. 2002; Dc Troia et al. 2003). 
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Peebles & Ratra 2003), with h as a function of f^o computed accordingly. The values of were chosen 
to be consistent with current standard nucleosynthesis constraints (fish 2 — 0.006, 0.014, or 0.022, see, e.g., 
Peebles & Ratra 2003). To render the problem tractable, Ce's were determined for the central values of 
t and rish 2 , and for the two combinations of these parameters which most perturb the C/s from those 
computed at the central values (i.e., for the smallest to we used the smallest Vlsh 2 , and for the largest to 
we used the largest Fish 2 ). Specific parameter values are given in columns (1) and (2) of Tables 1-3, and 
representative CMB anisotropy power spectra can be seen in Figs. 1 and 2. 

The differences in the low-£ shapes of the C^'s in Figs. 1 and 2 are a consequence of the interplay between 
the "usual" (fiducial CDM) Sachs- Wolfe term and the "integrated" Sachs- Wolfe term in the expression for 
the CMB spatial anisotropy (see, e.g., Hu & Dodclson 2002 for a review). The relative importance of these 
terms is determined by the values of Slo and a. (These terms are not easy to estimate analytically, and 
we can not trust the inaccurate XCDM approximation and must instead rely on numerical computations 
because these terms depend in complicated ways on the way in which w changes from its value during the 
matter dominated epoch to its value at zero rcdshift.) More precisely, in dark energy models the gravitational 
potential decays with time at late time, so at late time the gravitational potential and its time derivative are 
of opposite sign. Consequently there is some anticorrelation between the usual and integrated Sachs- Wolfe 
terms. This anticorrelation more significantly reduces the C/s at larger Qq (where the integrated term is 
not too large) and larger a (because the transition from matter dominance to dark energy dominance occurs 
earlier). 

The normalization of the theoretical CMB anisotropy spectra are determined by comparing them to 
the DMR data using the likelihood analysis method of Gorski (1994). In this paper we utilize the DMR 
four-year co-added 53 and 90 GHz sky maps in galactic coordinates. We do not consider the DMR ecliptic 
coordinates data here; the small shifts in the inferred normalization amplitudes due to the small differences 
in the galactic- and ecliptic-coordinate maps arc quantified in Gorski et al. (1998). The 53 and 90 GHz 
maps are coadded using inverse-noise-variance weights. The main changes relative to the analysis of Gorski 
ct al. (1998) are that we use the DMR maps in the HEALPix pixelisation, and excise pixels likely to be 
contaminated by bright foreground emission on and near the Galactic plane (Banday et al. 1997) using a 
cut rederived from the DIRBE 140 /xm data. The DIRBE map was reconstructed in the HEALPix format 
using the publicly available Calibrated Individual Observations (CIO) files together with the DIRBE Sky 
and Zodi Atlas (DSZA, both products being described in Hauser et al. 1998), allowing the construction of a 
foil sky map corrected for the zodiacal emission according to the model of Kelsall ct al. (1998). We consider 
C/s to multipole I < 30. See Gorski ct al. (1998) for a more detailed description of the method. 

As mentioned, we excise sky-map pixels on and near the Galactic plane where foreground emission 
dominates the CMB. Following Gorski et al. (1998, also see Gorski et al. 1995, 1996) we quantify the extent 
to which residual high-latitude Galactic emission can modify our results in two ways: by performing all 
computations both including and excluding the observed sky quadrupole, and by performing all computations 
with and without correcting for emission correlated with the DIRBE 140 (im sky map, itself fitted to the 
data for each model. 11 



11 See Gorski et al. (1998) for a more detailed discussion. More recent discussions of foreground emissions may be accessed 
through Mukherjee et al. (2002, 2003b), Banday ct al. (2003), and Bennett et al. (2003b). 
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3. Results 
3.1. Results of Qrms-ps fitting 

The results of the DMR likelihood analysis are summarized in Figs. 3-9 and Tables 1-3. 

We have computed DMR data likelihood functions £(Q r ms-PSi a ) using three sets of C/s computed 
for values of (fl B h 2 , t ) = (0.014, 13 Gyr), (0.006, 11 Gyr), and (0.022, 17 Gyr). We have also considered four 
"different" DMR data sets, either accounting for or ignoring the correction for faint high-latitude foreground 
Galactic emission, and either including or excluding the quadrupole moment from the analysis. 

Variations in f^o and a have the largest effect on the inferred value of Qrms-PS from the DMR data. 
Changes in fish 2 , to, faint highdatitude Galactic emission treatment, and quadrupole moment treatment, 
have much smaller effects, consistent with the findings of Gorski et al. (1998). These results determine what 
we have chosen to display in Figs. 3-9 and Tables 1-3. 

Two representative sets of likelihood functions L(Q rms _ps, ^o) are shown in Figs. 3 and 4. Figure 3 
shows those derived while accounting for the correction for faint highdatitude foreground Galactic emission, 
and including the quadrupole moment in the analysis. Figure 4 shows the likelihood functions derived 
while ignoring the faint highdatitude foreground Galactic emission correction, and excluding the quadrupole 
moment from the analysis. These data sets represent two different ways of dealing with effects of foreground 
emission on the inferred Qrms-ps normalization (with the former resulting in a typically 5 % lower <3 rms -ps 
value). These two data sets span the range of normalizations inferred from our analysis here. 12 The other two 
data sets we have used here, accounting for the faint highdatitude foreground emission while excluding the 
quadrupole and not accounting for the faint highdatitude foreground emission while including the quadrupole, 
result in Q rm s-ps normalizations that are typically 1-2 % above the smaller value or below the larger value 
of the first two data sets. 

Tables 1-3 give the Qrms-ps central values and 1 and 2 a ranges, computed from the appropriate 
posterior probability density distribution function assuming a uniform prior. These tables are computed for 
the three cases with tt B h 2 = 0.014 and t = 13 Gyr (Table 1), fl B h 2 = 0.006 and t = H Gyr (Table 2), and 
fish 2 = 0.022 and to = 17 Gyr (Table 3). We present results for four representative values of a (= 0, 2, 4, 
and 6), and a range of values of flo spanning the interval from 0.1 to 1. <3 rms -ps values are given only for the 
two extreme data sets whose likelihood functions are shown in Figs. 3 and 4: the faint highdatitude Galactic 
foreground emission corrected, quadrupole included case, and the Galactic foreground emission uncorrected 
quadrupole excluded case. 

Figure 5 shows the ridge lines of maximum likelihood Q rms _ps value, as a function of fig, for two 
different values of a, for all four data sets used in this paper. The upper panels of Fig. 7 show the effect the 
small differences between data sets have on the conditional (fixed Clo, a, Qsh 2 , and to) likelihood function 
for <3 rms -ps. 

Tables 1-3 also illustrate the shift in the inferred Qrms-ps normalization amplitudes due to changes 
in to and fl B h 2 . Figure 6 shows the effects that varying to and fl B h 2 have on the ridge lines of maximum 
likelihood Qrms-ps as a function of fl - The lower panels of Fig. 7 show the effects varying to and Q, B h 2 on 
the conditional (fixed SIq and a) likelihood functions for <2 rms -ps for a given data set. On the whole, for 



12 The ecliptic-frame sky maps result in slightly larger inferred Q r ms-PS normalization values than the galactic-frame sky 
maps used here, see Gorski et al. (1998). 



- 6 - 



the CMB anisotropy spectra considered here, shifts in to and Qsh 2 have only a small effect on the inferred 
normalization amplitude. 

Future WMAP data (see, e.g., Bennett et al. 2003a) should be able to significantly constrain the 
dark energy scalar field model. Here we use the DMR normalized dark energy model predictions to semi- 
quantitatively constrain the cosmological parameters of this model. Since we do not attempt to determine 
precise constraints on these parameters in this paper, we do not combine results from analyses of all DMR 
data sets, as Gorski et al. (1998) did for the open model (Gott 1982; Ratra & Peebles 1995), to determine 
the most robust estimate of the Qrms-ps normalization amplitudes. 

Figures 8 and 9 show marginal likelihood functions for some of the cosmological model parameter 
values and DMR data sets considered here. Figure 8 shows the likelihood function L(f2 ,aO derived by 
marginalizing i(fio, a, Qrms-ps) over Qrms-ps with a uniform prior. Figure 9 shows conditional (on slices of 
constant a or ilo) projections of the two-dimensional likelihood functions of Fig. 8. The effects of Galactic 
foreground correction and inclusion of the quadrupole are important. While the DMR data qualitatively 
indicate that lower a and higher Q values are mildly favored, these data by themselves do not provide 
significant quantitative constraints on cosmological model parameter values. 



The fractional energy-density perturbation power spectrum P(k) (eq. [1]) was evaluated from a numer- 
ical integration of the linear perturbation theory equations of motion. Table 4 lists the P(k) normalization 
amplitudes A (eq. [1]) for a range of i7 and a values. 13 The normalization amplitude is quite insensitive to 
the values of f2#/i 2 and to used. It is more sensitive to the choice of DMR data set, and Fig. 10 shows repre- 
sentative P(k)'s normalized to Q rms _ps obtained by averaging those from the foreground-emission-corrected 
quadrupole-included analysis and the foreground-emission-uncorrected quadrupole-excluded analysis. Fig- 
ure 10 also shows recent observational determinations of the galaxy power spectrum on linear scales. The 
predicted matter power spectrum P(k) <x (Qrms-ps) 2 so the error in the DMR normalization of P(k) is ap- 
proximately twice that in the DMR determination of Qrms-ps (which is discussed below). As always, when 
comparing model predictions to observational measurements it is important to account for this normaliza- 
tion uncertainty, rather than basing conclusions about model viability solely on the "central" normalization 
value. 

The mean square linear mass fluctuation averaged over a sphere of coordinate radius \ is 



Tables 5-7 list predictions for (<5M/M)(8/i _1 Mpc) or erg, the square root of the above expression evalu- 
ated at 8/i _1 Mpc, for various values of a and Q . In Tables 5-7 we list the mean of the central values 
obtained from the forcground-emission-corrected quadrupole-included analysis and the foreground-emission- 
uncorrected quadrupole-excluded analysis. 

Gorski et al. (1998) show that the maximal 1 a fractional Qrms-ps ( an d so SM/M) uncertainty from the 
DMR data ranges between about 10 % and 12 %, depending on model parameter values, and the maximal 2 



13 Numerical noise and other small effects contribute to the tiny artificial differences between the various Qq = 1 model 
predictions in this and other tables. 



3.2. Computation of Large-Scale Structure Statistics 
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a uncertainty ranges between about 16 % and 19 %. As we have not done a complete analysis of all possible 
DMR data sets here, we adopt 10 % and 20 % as our 1 and 2 a DMR normalization uncertainties for the dark 
energy scalar field models, for the purpose of comparing model predictions to observational measurements. 

4. Observational Constraints on DMR-Normalized Models 

Given the dependence of the DMR likelihoods on the quadrupole moment of the CMB anisotropy, the 
DMR data do not meaningfully exclude any part of the (Oo ; a, to, Sl B h 2 ) parameter space for the dark 
energy scalar field model considered here. In this section we combine current observational constraints on 
cosmological parameters (as summarized in Peebles & Ratra 2003; Chen & Ratra 2003b; Bennett et al. 2003a) 
with the DMR-normalized model predictions to place semi-quantitative constraints on the range of allowed 
model-parameter values. 

We focus on the parameter values and predictions of Tables 5-7. These are computed for parameter 
values n B h 2 = 0.014 and t = 13 Gyr (Table 5), fl B h 2 = 0.006 and t = 11 Gyr (Table 6), and fl B h 2 = 0.022 
and to = 17 Gyr (Table 7), which span the two standard deviation range of interest: 

11 Gyr < t < 17 Gyr, (6) 
0.006 < fl B h 2 < 0.022, (7) 

with central values of to = 13 Gyr and Vl B h 2 = 0.014 (as summarized in Peebles & Ratra 2003; see Carretta 
ct al. 2000, Krauss & Chaboyer 2001, and Chaboyer & Krauss 2002 for t , and Buries, Nollett, & Turner 
2001 and Cyburt, Fields, & Olive 2001 for n B h 2 ). The recent WMAP measurements (Bennett ct al. 2003a) 
lie in the center of the above to range, 13.3 Gyr < to < 14.1 Gyr (2 a), and near the upper end of the above 
£l B h 2 range 0.021 < Q B h 2 < 0.024 (2 a), in good agreement with standard Big Bang nucleosynthesis theory 
and primordial deuterium abundance observations. 

The first columns in Tables 5-7 list f^o- A recent analysis of pre- WMAP constraints on Qq from a com- 
bination of dynamical, baryon fraction, power spectrum, weak lensing, and cluster abundance measurements 
results in the two standard deviation range 

0.2 < il n £ 0.35 (8) 

(Chen & Ratra 2003b), in striking accord with the Bennett et al. (2003a) WMAP estimate of 0.19 < fl < 
0.35 (2 a). This constraint on f2 significantly narrows the range of viable model-parameter values. 

From an analysis of all available measurements of the Hubble parameter prior to mid 1999, Gott ct 
al. (2001) find at two standard deviations 

0.6 < h < 0.74, (9) 

with central value h = 0.67. This is in good agreement with the WMAP estimate 0.65 < h < 0.79 (2 a) 
with central value h = 0.71 (Bennett et al. 2003a). This constraint on h also significantly narrows the range 
of viable model parameters. 

For an estimate of the rms linear fractional mass fluctuations at 8/i _1 Mpc we use, at 2 cr, 



0.72 < a s < 1.16, 



(10) 
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with erg = 0.94 as central value (Chen & Ratra 2003b), with an additional 2 a DMR normalization uncertainty 
of 20 % added in quadrature. 

Focussing on the numerical values of Tables 5-7, and combining the preceding constraints, there remains 
a small range of parameter space that is observationally viable. In particular, models with £lo ~ 0.3 and a 
ranging from to ~ 3 with S!b/i 2 = 0.014 and t n — 13 Gyr, are viable, while larger values of fish 2 = 0.022 
and t = 17 Gyr narrow the range of a significantly, barely allowing a constant A model with a = 0. 
Dark energy models that deviate significantly from Einstein's cosmological constant are not favored by the 
data, i.e., the data favor a closer to 0. A model with a — 3, as is still allowed, requires new physics near 
a characteristic energy scale ~ 10 3 GeV (Peebles & Ratra 1988, 2003), which should soon be probed by 
accelerator-based elementary particle physics experiments. 

5. Conclusion 

We have normalized the inverse-power-law-potential scalar field dark energy model by comparing the 
predicted CMB anisotropy in this model to that measured by the DMR experiment. In addition to effects 
associated with the DMR data — see Gorski et al. (1998) for a more detailed discussion — the model 
normalization is sensitive to the values of fio and a (the exponent of the inverse-power-law potential) but 
almost independent of the values of £Ib and h. 

The DMR data alone can not be used to constrain Q and a over the ranges considered here. Current 
cosmographic observations, in conjunction with current large-scale structure observations compared to the 
predictions of the DMR-normalizcd dark energy scalar field model, do however provide quite restrictive 
constraints on model parameter values. This is largely because <j% depends sensitively on the values of a and 
f2o (see, e.g., Douspis et al. 2002), more so than the positions of the peaks in the CMB anisotropy spectrum. 

As discussed in the previous section, there is a small range of model parameter values where evolving 
dark energy density models remain an observationally-viable alternative to a constant cosmological constant 
model, however, the data do not favor large deviations from the constant A model. It would therefore be of 
interest to examine more closely the supergravity-inspired generalization of the inverse-power-law-potential 
scalar field dark energy model (Brax & Martin 2000; Brax et al. 2000), which has a scalar field potential 
that results in a slower evolution of the dark energy density, more like a constant A. 

It is possible that current or future WMAP (see, e.g., Bennett et al. 2003a) and other CMB anisotropy 
data will significantly circumscribe this option (see, e.g., Baccigalupi et al. 2002; Mukhcrjcc et al. 2003a; 
Caldwell & Doran 2003), although it should be noted that the constraints derived here are quite restrictive 
since as is very sensitive to the values of a and fig. Among other experiments, the proposed SNAP space 
mission to measure the redshift-magnitude relation out to a rcdshift of two should also provide interesting 
constraints on evolving dark energy (see, e.g., Podariu, Nugent, & Ratra 2001a; Wang & Lovelace 2001; 
Weller & Albrecht 2002; Gerke & Efstathiou 2002; Eriksson & Amanullah 2002; Rhodes et al. 2003). 

We acknowledge the advice and assistance of P. Brax, J. Lesgourgues, and J. Martin and thank the 
referee, M. Vogeley, for helpful comments. PM and BR acknowledge support from NSF CAREER grant AST- 
9875031. This work was partially carried out at the Jet Propulsion Laboratory of the California Institute of 
Technology, under a contract with the National Aeronautics and Space Administration. 
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Q921.28 22.76 
yz 18. 64 17.44 


20. 


1C21.56 23.12 
iD 18. 84 17.56 


21 


0922. 80 24.40 
• Jz 19.96 18.64 


21 


1R 22.68 24.28 
■ iD 19.76 18.40 


22. 


4023.96 25.64 
* u 20.96 19.56 


0.5 





.62 


17, 


/tnl8.60 19.88 
* u 16.28 15.24 


18. 


oel9.60 20.92 
• 5D 17.24 16.08 


17 


r 9 18.72 20.00 
• OZ 16.40 15.32 


18. 


40.19. 72 21.08 
*°17.32 16.20 


18. 


9019.44 20.80 
zu 17.00 15.88 


19 


9n20. 52 21.92 
■ zu 18.00 16.80 


18 


n.920. 20 21.64 
■ yz 17.68 16.52 


20. 


0021-36 22.80 
uu 18.72 17.48 


0.6 





.59 


17. 


i 9I8.28 19.52 
lz 16.04 15.00 


18. 


04 19.24 20.56 
u *16.92 15.84 


16 


4S.17.60 18.84 
■^°15.44 14.44 


17. 


/inl8.52 19.80 
* u 16.28 15.24 


16. 


rol7. 64 18.84 
oz 15.48 14.44 


17 


/inl8.56 19.84 
■* u 16.32 15.28 


16 


onl7.92 19.16 
■ ou 15.72 14.68 


17, 


«ol8.88 20.12 
D °16.60 15.52 


0.7 





.56 


16. 


QR18.12 19.36 
yo 15.88 14.84 


17. 


ool9.08 20.36 
°°16.76 15.68 


15 


nel7.00 18.20 
■ yo 14.92 13.96 


16. 


onl7.92 19.08 
•° u 15.76 14.76 


15. 


4416.48 17.60 
^14.48 13.52 


16 


9417.32 18.48 
■ z ^15.24 14.28 


15 


1 0I6.I6 17.24 
■ lo 14.20 13.28 


15, 


0916.96 18.12 
yz 14.96 14.00 


0.8 





.54 


16. 


QO18.08 19.32 
yo 15.84 14.80 


17. 


ool9.04 20.32 
•°°16.72 15.68 


15 


ool6.96 18.08 
■°°14.88 13.88 


16. 


7917.80 19.00 
' z 15.68 14.68 


15. 


O»16.08 17.16 
u °14.12 13.20 


15 


00I6.88 18.04 
'°°14.88 13.96 


14 


4ol5.4416.48 
■^°13.60 12.72 


15. 


94 16.24 17.32 
z ^14.32 13.40 


0.9 





.52 


16. 


QR18.12 19.36 
yo 15.88 14.84 


17. 


0919.08 20.36 
yz 16.76 15.72 


16 


OO17.40 18.60 
■ oz 15.28 14.28 


17. 


1 fi 18.32 19.56 
i0 16.12 15.08 


15. 


a a 16.72 17.84 
°*14.68 13.72 


16 


4»17.56 18.72 
■* 15.48 14.48 


15 


94l6.2417.32 
■ z *14.28 13.32 


16. 


0017.08 18.20 
uu 15.04 14.08 


1 





.50 


17. 


n/| 18.20 19.48 
u *15.96 14.92 


18. 


nnl9-20 20.48 
uu 16.88 15.76 


17 


n ol8. 24 19.52 
■ uo 16.00 14.96 


18. 


n/|19-20 20.52 
u *16.92 15.80 


17. 


n/| 18.20 19.44 
uq 45.96 14.88 


17 


qg19. 16 20.44 
■ yD 16.84 15.76 


17 


04 18.20 19.48 
• uq 45.96 14.92 


18. 


0019 20 20.48 
uu 16.88 15.76 



a Thc tabulated Qrms-PS values are determined from the conditional likelihood function at fixed fio and a. At each fi and a, the first of the five entries in 
each of columns (3)-(10) is the maximum likelihood value, the first (vertical) pair delimits the 68.3 % (1 a) highest posterior density range, and the second 
(vertical) pair delimits the 95.5 % (2 a) highest posterior density range. 

b Accounting for (Yes), or ignoring (No) the correction for faint high-latitude foreground Galactic emission. 



Table 2. Q rms _ PS Values for the t = 11 Gyr, n B h 2 = 0.006 Models a 



a=0 a=2 a=4 a=6 



G.C. b : Yes No Yes No Yes No Yes No 

^min* 2 3 2 3 23 23 



il 




h 




Qrms-PS 

( M K) 




Qrms— PS 

(jjK) 




Qrms-PS 

( M K) 




Qrms-PS 




Qrms-PS 

(/xK) 




Qrms-PS 

( M K) 




Qrms — PS 

( M K) 




Qrms-PS 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




(8) 




(9) 




(10) 


0.1 


1 


.13 


22. 


7024.36 26.12 
' z 21.16 19.68 


24, 


1925.80 27.64 
iz 22.52 20.96 


24 


QR26.84 28.84 
■ yD 23. 24 21.56 


25. 


QR28.52 30.60 
yD 24.76 23.04 


26 


qr28. 96 31.16 
■ yo 25.04 23.24 


28. 


7930.80 32.96 
' z 26.76 24.84 


28 


7R30.92 33.16 
■'°26.76 24.80 


30 


fi/|32.76 35.12 
°*28.56 26.52 


0.2 





.96 


20. 


oo21.80 23.36 
'°°19.00 17.72 


21, 


CR23.04 24.68 
do 20.16 18.84 


22 


R024.32 26.12 
•°°21.12 19.64 


24. 


ns25. 80 27.64 
uo 22.48 20.92 


24 


CR26.36 28.32 
' do 22. 84 21.24 


26. 


1927.96 30.00 
lz 24.32 22.64 


26 


no27.96 30.04 
' uu 24.20 22.48 


27. 


co29. 68 31.84 
°°25.80 24.04 


0.3 





.86 


18. 


QR 20. 28 21.72 
yD 17.72 16.52 


20. 


n/ |21.40 22.88 
u ^18.76 17.56 


20 


7C22.24 23.84 
' ' °19.36 18.04 


22, 


nn23. 52 25.16 
uu 20.56 19.16 


22 


a a 24.04 25.80 
■^20.88 19.44 


23 


on25.48 27.28 
° u 22.24 20.68 


23 


a a 25.36 27.24 
'°*22.04 20.48 


25. 


1926.88 28.84 
iz 23. 48 21.84 


0.4 





.79 


18. 


1919.36 20.72 
iz 16.92 15.80 


19, 


1 920.44 21.84 
1Z 17.92 16.76 


19 


n/l 20.36 21.80 
■ u *17.76 16.56 


20, 


1921.52 23.04 
iz 18.84 17.60 


20 


o R 21.80 23.32 
■^°19.00 17.68 


21, 


r 9 23.04 24.64 
OZ 20.16 18.80 


21 


o R 22.88 24.52 
■•J D 19.96 18.56 


22, 


ca 24.20 25.92 
D *21.16 19.76 


0.5 





.74 


17, 


Rn 18.80 20.12 
DU 16.48 15.36 


18. 


c R 19.80 21.16 
OD 17.40 16.24 


17 


00I8.88 20.24 
• D °16.52 15.44 


18. 


R ol9. 92 21.28 
D5 17.48 16.32 


18 


OR19.64 21.00 
•' 5D 17.16 16.04 


19. 


/in20.72 22.16 
* u 18.16 16.96 


19 


1 920.44 21.88 
■ lz 17.88 16.64 


20. 


9021.60 23.08 
zu 18.92 17.68 


0.6 





.70 


17, 


r>ol8.44 19.72 
z °16.16 15.08 


18. 


9/119.44 20.76 
z ^17.08 15.96 


16 


en!7.76 19.00 
■ DU 15.56 14.52 


17. 


C9I8.72 20.00 
oz 16.44 15.36 


16 


ool7.80 19.04 
■ DC, 15.60 14.56 


17, 


cnl8 .76 20.04 
DU 16.48 15.40 


16 


QR18.12 19.36 
• yD 15.88 14.84 


17, 


ool9.08 20.36 
°°16.76 15.68 


0.7 





.66 


17, 


19I8. 24 19.56 
lz 16.00 14.96 


18, 


n4 19-24 20.56 
u ^16.92 15.80 


16 


n«17.20 18.36 
' u °15.08 14.04 


16, 


QR18.08 19.32 
yo 15.88 14.88 


15 


Rn!6.64 17.76 
■ OU 14.60 13.64 


16. 


/inl7.48 18.68 
^"15.36 14.40 


15 


90I6.32 17.44 
■ zo 14.32 13.40 


16, 


nol7.16 18.28 
u °15.12 14.12 


0.8 





.64 


17, 


n«18. 24 19.48 
uo 16.00 14.92 


18, 


nnl9-20 20.52 
uu 16.88 15.80 


16 


0/117.12 18.28 
u *15.00 14.00 


16, 


00I8.OO 19.20 
°°15.84 14.80 


15 


9/116.2417.36 
• z ^14.24 13.32 


16 


nn!7.08 18.20 
uu 15. 04 14.04 


14 


a a 15.60 16.68 
•°^13.72 12.84 


15. 


/inl6.40 17.48 
* u 14. 44 13.52 


0.9 





.61 


17. 


1 9I8.32 19.56 
lz 16.04 15.00 


18. 


n ol9.28 20.60 
uo 16.92 15.84 


16 


/|/|17.56 18.76 
•^15.40 14.36 


17. 


09I8.48 19.72 
oz 16.24 15.20 


15 


on!6.88 18.00 
■ OU 14.80 13.80 


16. 


a a 17.72 18.92 
'°*15.60 14.60 


15 


o«16.40 17.48 
■ OD 14.40 13.44 


16. 


1 fi 17.20 18.40 
i0 15.16 14.20 


1 





.59 


17. 


9018.40 19.64 
zu 16.12 15.04 


18. 


1R19.36 20.68 
iD 17.00 15.92 


17 


9/118.40 19.68 
■ zq 46.12 15.04 


18, 


9O19-40 20.72 
zu 17.04 15.96 


17 


1 R18.36 19.60 
■ iD 16.08 15.00 


18, 


1919.32 20.64 
iz 16.96 15.88 


17 


9O18.40 19.64 
■ zu 16.12 15.04 


18, 


1R19.36 20.68 
iD 17.00 15.92 



a Thc tabulated Qrms-PS values are determined from the conditional likelihood function at fixed fio and a. At each fi and a, the first of the five entries in 
each of columns (3)-(10) is the maximum likelihood value, the first (vertical) pair delimits the 68.3 % (1 a) highest posterior density range, and the second 
(vertical) pair delimits the 95.5 % (2 a) highest posterior density range. 

b Accounting for (Yes), or ignoring (No) the correction for faint high-latitude foreground Galactic emission. 



Table 3. Q rms _ PS Values for the t = 17 Gyr, n B h 2 = 0.022 Models a 



a=0 a=2 a=4 a=6 



G.C. b : Yes No Yes No Yes No Yes No 

^min* 2 3 2 3 23 23 







h 




Qrms-PS 

( M K) 




Qrms— PS 

(jjK) 




Qrms-PS 




Qrms-PS 

( M K) 




Qrms-PS 




Qrms-PS 

( M K) 




Qrms — PS 

( P K) 




Qrms-PS 


(1) 


(2) 




(3) 




(4) 




(5) 




(6) 




(7) 




(8) 




(9) 




(10) 


0.1 





.73 


22. 


9^23.76 25.44 
zu 20.72 19.28 


23 


ro25. 12 26.88 
oz 21.96 20.52 


24 


,1026.28 28.20 
'* o 22.80 21.20 


26. 


no27.88 29.84 
uu 24.28 22.60 


26 


9,128.20 30.28 
■ z *24. 44 22.72 


27 


0929.92 32.04 
■ yz 26.04 24.28 


27. 


oa 29. 92 32.08 
°*25.92 24.08 


29. 


R031.72 33.88 
°°27. 68 25.76 


0.2 





.62 


19. 


oo21.24 22.76 
°°18.56 17.32 


21. 


n n22.44 23.96 
uu 19.68 18.40 


22 


9023.84 25.60 
' zo 20.76 19.32 


23 


Rn 2 5. 24 27.04 
OU 22.04 20.56 


24 


n/| 25. 08 27.68 
■ u ^22.40 20.84 


25 


ce27.36 29.28 
• do 23. 84 22.24 


25. 


OR27.20 29.20 
'°°23.60 21.96 


26. 


QR28.88 30.92 
yo 25.16 23.44 


0.3 





.55 


18. 


ar\W.88 21.24 
OU 17.40 16.24 


19. 


a a 20.96 22.40 
°%8.40 17.20 


20 


QR21.76 23.36 
■ OD 19.00 17.72 


21. 


CR23.04 24.60 
OD 20.16 18.80 


21 


QR23.52 25.24 
■ yo 20.48 19.08 


23 


9024.92 26.64 
• zo 21.76 20.28 


23 


n o24.76 26.52 
uo 21.56 20.04 


24. 


C926.20 28.08 
oz 22. 88 21.36 


0.4 





.51 


17. 


7R19.00 20.28 
' D 16. 64 15.52 


18. 


7R20.00 21.36 
•' D 17.56 16.44 


18 


CO20.00 21.36 
■ D5 17.48 16.28 


19 


7R21.08 22.52 
' D 18.48 17.28 


19 


Q921.32 22.84 
■ yz 18.64 17.36 


21 


0O22.52 24.08 
■ uo 19.72 18.44 


20. 


oo22.36 23.92 
°°19.52 18.16 


22. 


1923. 60 25.24 
■ iz 20.68 19.28 


0.5 





.48 


17. 


o n 18.36 19.64 
zu 16.12 15.04 


18. 


1R 19.36 20.64 
1D 17. 04 15.92 


17 


ool8.48 19.76 
•' 5Z 16.20 15.12 


18. 


9/119.48 20.80 
z4 17.12 16.00 


18 


nnl9-24 20.56 
■ uu 16. 84 15.72 


19 


nf|20.28 21.64 
■ uu 17.80 16.64 


18. 


eol9.92 21.32 
• D0 17.48 16.28 


19. 


7921.04 22.44 
■ ' z 18.48 17.28 


0.6 





.45 


16. 


QR18.12 19.36 
yo 15.92 14.88 


17. 


ool9.08 20.36 
""le.SO 15.72 


16 


ool7.40 18.60 
■ oz 15.28 14.28 


17. 


9nl8.32 19.56 
zu 16.12 15.08 


16 


o R 17. 44 18.68 
■°°15.32 14.32 


17 


9/118.36 19.60 
• z %6.16 15.12 


16. 


rccl7.68 18.92 
OD 15.52 14.52 


17. 


40I8.6O 19.84 
*°16.36 15.32 


0.7 





.43 


16. 


7R17.92 19.12 
• ' °15.72 14.68 


17. 


ro18.84 20.12 
°°16.56 15.52 


15 


onl6.88 18.04 
• ou 14.84 13.84 


16. 


a a 17.76 18.92 
•°^15.64 14.60 


15 


90I6.28 17.36 
■ zo 14.32 13.36 


16 


A/1 17.08 18.24 
• u ^15.04 14.12 


14. 


QR15.92 17.00 
yo 14.00 13.12 


15. 


7916.72 17.84 
' ' z 14.76 13.80 


0.8 





.41 


16. 


onl7.92 19.12 
ou 15.72 14.68 


17. 


fiol8.84 20.12 
°°16.60 15.52 


15 


79I6.8O 17.92 
■ ' z 14.76 13.76 


16. 


Cfil7. 64 18.84 
,JU 15.52 14.56 


14 


QR15.96 17.04 
■ yo 14.04 13.12 


15 


7916.76 17.88 
•' z 14.76 13.84 


14. 


09I5.28 16.28 
oz 13.44 12.60 


15. 


n«16. 04 17.08 
uo 14.16 13.24 


0.9 





.40 


16. 


oa 17.96 19.20 
' *15.76 14.72 


17. 


7R18.92 20.16 
'°16. 64 15.56 


16 


iel7. 24 18.44 
■ iD 15.16 14.16 


17. 


nnl8.12 19.36 
uu 15.96 14.96 


15 


C9I6.56 17.68 
■ oz 14.56 13.60 


16 


o fi 17.40 18.56 
■ OD 15.32 14.36 


15. 


n/f 16.00 17.08 
u *14.08 13.16 


15. 


oa16. 84 17.92 
ou 14.84 13.88 


1 





.38 


16. 


oa 18.04 19.28 
■ oq: 15.84 14.80 


17. 


Q/i 19.00 20.28 
°*16.72 15.64 


16 


Q0I8.O8 19.32 
■ yz 15.88 14.84 


17. 


O/i 19.00 20.28 
■°*16.72 15.64 


16 


Q9I8. 04 19.28 
■ yz 15.84 14.80 


17 


on!9-00 20.28 
■ ou 16.72 15.64 


16. 


Q9I8.O8 19.32 
• yz 15.88 14.84 


17. 


oa 19.04 20.32 
°*16.76 15.68 



a Thc tabulated Qrms-PS values are determined from the conditional likelihood function at fixed fio and a. At each fi and a, the first of the five entries in 
each of columns (3)-(10) is the maximum likelihood value, the first (vertical) pair delimits the 68.3 % (1 a) highest posterior density range, and the second 
(vertical) pair delimits the 95.5 % (2 a) highest posterior density range. 

b Accounting for (Yes), or ignoring (No) the correction for faint high-latitude foreground Galactic emission. 
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Tabic 4. Fractional Energy-Density Perturbation Power Spectrum Normalization Factor Ah' 





a=0 


a=2 


a=A 


a=6 




Ah 4 


Ah 4 


Ah 4 


Ah 4 




(10 5 Mpc 4 ) 


(10 5 Mpc 4 ) 


(10 5 Mpc 4 ) 


(10 5 Mpc 4 ) 


(i) 


(2) 


(3) 


(4) 


(5) 


0.1 


153 (28.4) 


307 (47.0) 


525 (69.1) 


802 (92.9) 


0.2 


62.9 (14.6) 


125 (23.2) 


208 (33.1) 


303 (42.9) 


0.3 


36.7 (9.82) 


63.2 (14.1) 


101 (19.3) 


141 (24.2) 


0.4 


24.8 (7.28) 


34.3 (9.14) 


50.3 (11.7) 


66.8 (14.1) 


0.5 


18.0 (5.64) 


19.9 (6.14) 


25.0 (7.16) 


30.8 (8.13) 


0.6 


13.8 (4.48) 


12.4 (4.31) 


12.9 (4.49) 


13.9 (4.69) 


0.7 


11.0 (3.62) 


8.62 (3.21) 


7.56 (3.01) 


6.92 (2.87) 


0.8 


8.99 (2.96) 


6.89 (2.59) 


5.51 (2.30) 


4.55 (2.06) 


0.9 


7.46 (2.45) 


6.29 (2.24) 


5.28 (2.05) 


4.58 (1.88) 


1 


6.29 (2.05) 


6.31 (2.05) 


6.27 (2.05) 


6.29 (2.05) 



"■Normalized to the mean of the central Qrms-ps in the foreground- 
emission-corrcctcd quadrupole-included and the foreground-emission-un- 
corrected quadrupole-excluded cases (normalized to Q rm s-ps = 10 /-tK); 
scale as (Qrms-ps) 2 - These are for the t = 13 Gyr and fish 2 = 0.014 
models and are insensitive to the values of these parameters. 



Table 5. (8M/M)(8h~ 1 Mpc) for the t = 13 Gyr, n B h 2 = 0.014 Models a 







a=0 


a=2 




a=6 




h 






CT8 


0-8 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0.1 


0.96 


0.632 


0.295 


0.169 


0.107 


0.2 


0.81 


0.888 


0.509 


0.332 


0.234 


0.3 


0.72 


1.02 


0.670 


0.479 


0.363 


0.4 


0.67 


1.10 


0.798 


0.617 


0.495 


0.5 


0.62 


1.14 


0.906 


0.745 


0.633 


0.6 


0.59 


1.17 


0.991 


0.863 


0.771 


0.7 


0.56 


1.18 


1.06 


0.972 


0.905 


0.8 


0.54 


1.18 


1.11 


1.06 


1.02 


0.9 


0.52 


1.18 


1.15 


1.12 


1.11 


1 


0.50 


1.17 


1.17 


1.17 


1.17 


a The mean 


of the central 


values 


obtained from 


the fore- 



ground-emission-corrected quadrupole-included analysis and the 
foreground-emission-uncorrected quadrupole-excluded analysis 



Table 6. (SM/M)(8h- 1 Mpc) for the t = 11 Gyr, VL B h 2 = 0.006 Models a 







a=0 


a=2 




a~6 


^0 


h 




CT8 


CT8 


CT8 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0.1 


1.13 


0.964 


0.477 


0.283 


0.185 


0.2 


0.96 


1.26 


0.746 


0.498 


0.356 


0.3 


0.86 


1.41 


0.945 


0.685 


0.522 


0.4 


0.79 


1.49 


1.10 


0.854 


0.689 


0.5 


0.74 


1.52 


1.22 


1.01 


0.862 


0.6 


0.70 


1.54 


1.31 


1.15 


1.03 


0.7 


0.66 


1.54 


1.39 


1.28 


1.19 


0.8 


0.64 


1.53 


1.44 


1.38 


1.33 


0.9 


0.61 


1.51 


1.47 


1.44 


1.44 


1 


0.59 


1.49 


1.49 


1.49 


1.49 


a The mean 


of the central 


values 


obtained from 


the fore- 



ground-cmission-corrected quadrupole-included analysis and the 
foreground-emission-uncorrected quadrupole-excluded analysis. 
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Table 7. (6M/M){8h,- 1 Mpc) for the t = 17 Gyr, n B h 2 = 0.022 Models a 







a=0 


a=2 


a=A 


a=6 




h 






CT8 


CT8 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0.1 


0.73 


0.273 


0.107 


0.0540 


0.0316 


0.2 


0.62 


0.453 


0.241 


0.149 


0.101 


0.3 


0.55 


0.562 


0.352 


0.244 


0.179 


0.4 


0.51 


0.631 


0.445 


0.336 


0.265 


0.5 


0.48 


0.675 


0.523 


0.425 
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a The mean of the central values obtained from the foreground- 
emission-corrected quadrupole-included analysis and the foreground- 
emission-uncorrected quadrupolc-cxcludcd analysis. 
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Fig. 1. — CMB anisotropy multipolc coefficients. Panels a) and d) in the first row, b) and e) in the second 
row, and c) and /) in the third row correspond to models with a = 0, 4, and 8, respectively. Panels a), 6), 
and c) in the left column show multipoles for a model with to = 13 Gyr and Slsh 2 = 0.014. The coefficients 
are normalized relative to the Cg amplitude, and different values of fio are offset from each other to aid 
visualization. In each of these panels, at I = 9, the 11 curves, in ascending order, correspond to models 
with = 0.05, 0.1, 0.2, . . . , 0.9, 1.0. The three panels in the right column show models with f2o = 0.2 (solid 
curves) and 0.4 (dashed curves). Each set of three curves in these panels correspond, in ascending order on 
the right vertical axis of each panel, to models with parameter values (to, 0.sh?) = (11 Gyr, 0.006), (13 Gyr, 
0.014), and (17 Gyr, 0.022), respectively. 
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Fig. 2. — CMB anisotropy multipolc coefficients. Panels a) and d) in the first row, b) and e) in the second 
row, and c) and /) in the third row correspond to models with Oo = 0.05, 0.3, and 0.6, respectively. Panels 
a), b), and c) in the left column show multipoles for a model with to = 13 Gyr and Slsft 2 = 0.014. The 
coefficients are normalized relative to the Cg amplitude, and different values of a are offset from each other 
to aid visualization. In each of these panels, at i = 9, the 9 curves, in ascending order, correspond to models 
with a = 0, 1, 2, . . . , 7, 8. The three panels in the right column show models with a = 4 (solid curves) and 1 
(dashed curves). Each set of three curves in these panels correspond, in ascending order on the right vertical 
axis of each panel, to models with parameter values (to, ^Ish 2 ) = (11 Gyr, 0.006), (13 Gyr, 0.014), and (17 
Gyr, 0.022), respectively. 



Fig. 3. — Likelihood functions L(<3rms-ps, ^o) (arbitrarily normalized to unity at the highest peak) derived 
from a simultaneous analysis of the DMR 53 and 90 GHz galactic-frame data, corrected for faint high-latitude 
foreground Galactic emission, and including the quadrupole moment in the analysis. These are for the model 
with to = 13 Gyr and fish 2 = 0.014. Panel a) is for a — and panel b) for a = 4. 




Fig. 4. — Likelihood functions L(Q rms „ps, Sl ) (arbitrarily normalized to unity at the highest peak) derived 
from a simultaneous analysis of the DMR 53 and 90 GHz galactic-frame data, ignoring the correction for 
faint high-latitude foreground Galactic emission, and excluding the quadrupole moment from the analysis. 
These are for the model with to = 13 Gyr and 0_b/i 2 = 0.014. Panel a) is for a = and panel b) for a = 4. 
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Fig. 5. — Ridge lines of the maximum likelihood Qrms-ps value as a function of fi , for a) a = and b) 
a = 4. These are for the model with to = 13 Gyr and Vtsh 2 = 0.014. Four curves are shown in each panel: 
solid (dashed) curves correspond to accounting for (ignoring) the faint high-latitude foreground Galactic 
emission correction, while heavy (light) curves include (exclude) the quadrupolc moment in the analysis. 




Fig. 6. — Ridge lines of the maximum likelihood Q rm s-ps value as a function of f^o, for a = 2. Six curves are 
shown. These correspond to parameter values (to, £lsh 2 ) = (11 Gyr, 0.006) (dotted curves), (13 Gyr, 0.014) 
(solid curves), and (17 Gyr, 0.022) (dashed curves). Heavy curves are determined by accounting for faint 
high-latitude foreground Galactic emission and including the quadrupole moment in the analysis, while light 
curves ignore the correction for faint high-latitude foreground Galactic emission and exclude the quadrupolc 
moment from the analysis. 



- 23 - 




15 20 25 30 15 20 25 30 

Q rms _ PS (/xK) Q rm5 _p S (,uK) 



Fig. 7. — Conditional likelihood densities for Q rms -ps, derived from _L((3 rms _ps, Oo, a) (which are normalized 
to be unity at the peak, for each DMR data set and set of model-parameter values). Left column panels 
a) and b) are for Q = 0.2 and right column panels c) and d) correspond to £1 = 0.4. Six curves are 
shown in each panel. Solid, dotted, and dashed curves correspond to a = 0, 2, and 4, respectively. In each 
panel, heavy curves are for a nominal model with to — 13 Gyr and f^ft 2 — 0.014 and are derived from 
the foreground-corrected quadrupole-includcd analysis. In the two upper panels the three light curves differ 
from the nominal heavy curves by not accounting for the faint high-latitude foreground Galactic emission 
correction and by excluding the quadrupolc from the analysis. In the two lower panels the three light curves 
differ from the nominal ones by corresponding to models with i = H Gyr and Os/i 2 = 0.006 (panel b) on 
the left) and with t = 17 Gyr and il B h 2 = 0.022 (panel d) on the right). 
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Fig. 8. — Likelihood functions L(a, CIq) (arbitrarily normalized to unity at the highest peak) derived by 
marginalizing L(Q rms -ps, a,flo) over Qrms-ps with a uniform prior. These are for the model with to = 
13 Gyr and fish 2 = 0.014. Panel a) ignores the correction for faint high-latitude foreground Galactic 
emission and excludes the quadrupole moment from the analysis, while panel b) corrects for faint high- 
latitude foreground Galactic emission and includes the quadrupole moment in the analysis. 
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Fig. 9. — Conditional likelihood functions derived from L(a, fio) (arbitrarily normalized to unity at the 
highest peak), which is derived by marginalizing L(Q r ms-PS, a , ^o) over <3 rms -ps with a uniform prior. 
These are for the model with t = 13 Gyr and fish 2 = 0.014. Panels a) and c) in the left column ignore the 
correction for faint high-latitude foreground Galactic emission and exclude the quadrupole moment from the 
analysis. Panels b) and d) in the right column account for the correction for faint high-latitude foreground 
Galactic emission and include the quadrupole moment in the analysis. The top two panels show L(Q ) on 
four constant a slices at a = (solid curves), 2 (dotted curves), 4 (dashed curves), and 8 (dot-dashed curves). 
The bottom two panels show L(a) on six constant S! slices at Clo = 1 (solid curves), 0.5 (dotted curves), 0.4 
(short dashed curves), 0.3 (dot-dashed curves), 0.2 (dash-three-dotted curves), and 0.1 (long dashed curves). 
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Fig. 10. — Fractional energy-density perturbation power spectra P(k) as a lunction of wavenumber k. These 
are normalized to the mean of the central Qrms-PS values from the foreground-corrected quadrupolc- included 
analysis and the foreground-uncorrected quadrupole-excluded analysis (for the corresponding values of the 
cosmological parameters). Panels a), b), and c) correspond to models with parameter values (to, fish 2 ) = 
(13 Gyr, 0.014), (11 Gyr, 0.006), and (17 Gyr, 0.022), respectively. Each panel shows 12 power spectra at 
three values of a = (solid curves), 2 (dotted curves), and 4 (dashed curves) for four values of Oo = 1, 0.4, 
0.2, and 0.1 in ascending order on the left axis of each panel. The triangles represent the linear real space 
power spectrum of the IRAS Point Source Catalogue Redshift Survey (PSCz) as estimated by Hamilton, 
Tcgmark, & Padmanabhan (2000). Upper limits are not plotted. The squares represent the redshift-space 
galaxy power spectrum estimated from the 2dFGRS optical galaxy data (Percival et al. 2001), again on 
linear scales, comparison with the shapes of curves being valid in the case of linear bias. 



